Generation of entanglement in regular systems 
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We study dynamical generation of entanglement in bipartite quantum systems, characterized 
by purity (or linear entropy), and caused by the coupling between the two subsystems. Explicit 
semiclassical theory of purity decay is derived for integrable classical dynamics of the uncoupled 
system, and for localized (general Gaussian wave-packet) initial states. Purity decays as an algebraic 
function of time x strength of perturbation, independently of the Planck's constant. 



PACS numbers: 3.65. Ud, 3.67.Mn, 3.65.Sq 
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Detailed understanding of entanglement, being one of 
the most distinct features of the quantum world, is an 
issue of high importance, particularly in view of recent 
efforts to build quantum devices that will manipulate 
(pure states of) individual quantum systems. The loss 
of control over the entanglement, e.g. decoherence, in 
such a device is one of the major obstacles that we have 
to overcome. 

In the present paper we are going to study dynami- 
cal generation of entanglement in bipartite systems. Ini- 
tially separable pure state will get entangled due to the 
coupling between two subsystems. Here we consider sys- 
tems where the uncoupled part of the Hamiltonian in 
both subsystems generates regular (integrable) dynamics 
in the classical limit. The motivation to study entan- 
glement generation in systems with a regular uncoupled 
dynamics comes from the fact that such systems are quite 
common both in experiments and as theoretical models. 
For instance, if the uncoupled system consists of a num- 
ber of uncoupled one degree of freedom (DOF) systems 
then it is integrable. Such is the case in various proposals 
for quantum computation, e.g. ion traps. Further, ex- 
perimentally realizable Jaynes-Cummings model, where 
decoherence for cat states [lj has actually been exper- 
imentally measured, is also an integrable system. Still 
further, a standard model of decoherence Q consists of 
an infinite number of harmonic oscillators. If the bath 
consists of a finite number of harmonic oscillators this 
falls under the domain of our theory. Recently |3( it has 
been pointed out that the decoherence for truly macro- 
scopic superposition is so fast that the usual master equa- 
tion approach is not valid anymore. On this very short 
"instantaneous" time scale any system will effectively be- 
have as a regular one (i.e. correlations do not decay yet). 
The semiclassical decay of purity has recently been dis- 
cussed by Jacquod Q , however his results in the regular 
domain do not agree with our findings. See Q for some 
other related recent results. 

Time evolution of the system will be governed by a 
Hamiltonian 

H = H + 5-V, H = H A ® 1 B + 1 A ®H B , (1) 

where Hq is an uncoupled part of the Hamiltonian and 
V is the coupling between the two subsystems respon- 
sible for the generation of entanglement. The strength 



of this coupling is given by a dimensionless parameter 
5. We will use subscripts "A" and "B" to denote two 
subsystems. The state of the whole system at time t is 
simply \ijj{t)) = U(t)\y^(0)) , with a unitary propagator 
U(t) = exp (— \Ht/K). Let us define time-averaged cou- 
pling 



V 



lira — 



dtV(t), 



(2) 



where V(t) is the coupling operator in the interaction 
picture, V(t) = U^(t)VU Q (t), U (t) = exp (-iH t/K), i.e. 
propagated with an uncoupled part of the Hamiltonian. 
We shall assume a situation, typical for a regular Hq, 
where V is a non-trivial operator, different from zero or 
a multiple of identity @ . We wish to stress that pertur- 
bation V will typically break the integrability of H, and 
that our results reported below are not limited to values 
of small S. 

The entanglement between the two subsystems, for a 
pure state \ip(t)}, is characterized by a purity 



I(t)= tv A [p 2 A (t)}, PA (t) = tv B \p(t)] 



(3) 



where p(t) := \tp(t))(t/j(t)\. Iff purity I(t) is less than 1, 
the two subsystems are entangled, otherwise they are in a 
separable (product) state. Our initial state will always be 
a product one, \ip(0)) = |^a(0)> <8> \i/>b(0)), hence 1(0) = 
1. The initial states |Va,b(0)) will be Gaussian wave 
packets. Time dependence of purity I(t) will then tell 
us how fast the two subsystems get entangled due to the 
coupling V. 

Let us proceed with the calculation of purity decay 
I(t). We should observe that propagating the state back- 
wards in time with a separable - uncoupled dynamics Uq (t) 
does not change the value of purity, so p(t) in Eq. can 
be replaced by 

p M (t) = M{t)p{0)M\t), M(t) = ul{t)U{t), (4) 

where M(t) is the echo operator used in the theory of 
fidelity decay 0, 0- The matrix p M (t) represents the 



evolution of our pure state in the interaction picture, 
just explained above, the purity © is equal to 



As 



I(t) = tr A [{pf(t)y] , pf(t) = tr B [p M (i)] . (5) 
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An advantage of the representation (J5J over © is the 
fact that the echo operator M${t) is, unlike the forward 
evolution U (t), close to an identity for small 8 so one may 
use perturbative, or asymptotic expansions in S. We fol- 
low the approach of Ref. |7f and use the Baker-Campbell- 
Hausdorff formula e sv e sw = exp(S(V+W) + ^S 2 [V, W} + 
. . .) for continuous products, see e.g. |j, to simplify the 
expression for the echo operator M(t). The lowest order 
term in the exponential is jrfdtV(t). For times larger 
than some classical averaging time i avc , on which V @ 
converges, this term can be rewritten as j-Vt. The second 
order term in 5 can be shown to grow with time no faster 
than 5 2 t/h, and by induction higher orders can be esti- 
mated to grow as ~ S r t r ~ 1 /h. Therefore, provided only 
8 <C 1, higher orders in S can be neglected and we end 
up with a very simple expression for the echo operator 



M(t) 



_ e -i6tV/h 



(0) 



So the echo operator can be interpreted as the prop- 
agator with an effective Hamiltonian SV. We proceed 
with a semiclassical evaluation of the purity, a procedure 
completely analogous to a similar calculation for the fi- 
delity |10|. We use the notation in which small Latin 
letters denote classical limiting observables (e.g. Weyl 
symbols) of the corresponding operators denoted by cap- 
ital Latin letters. For example, let j = (jaJb) denote 
a d = c£a + <fe dimensional vector of classical canonical 
actions of the completely integrable uncoupled classical 
Hamiltonian ho = h\ + h^. dj± and c?b are the number 
of DOF of the subsystems "A" and "B" , respectively. In 
quantum mechanics, one has a vector of mutually com- 
muting action operators J, with a common set of eigen- 
vectors, denoted by a multi-index n £ Z d of quantum 
numbers: J|n) = fi(n + a)|n) w /m|n) where a are the 
Maslov indices. Here and below means equal in the 
leading order in h. The purity (JSJ can now be written as 
a sum over d— dimensional lattice of quantum numbers, 
using the fact since V commutes with Hq it is diagonal 
in the basis |n), and in the leading semiclassical order (in 
h) wc can replace the summation by an integral over the 
classical action space. Further, we replace the operator 
V by its classical limit v(j), which is a conserved quan- 
tity so it is a function of d classical actions j only. Let us 
denote by p(j) = PA(jA)PB(jB) the classical limit of the 
initial density (n|p(0)|n). For our initial product state of 
two wave-packets each of the two densities is a Gaussian 



Pa(ja) = Cexp{-Q*-jt)K& -£)/*}, (7) 

where a subscript "a" takes values "A" or "B", de- 
pending on the subsystem, j* is the position of the 
initial packet, A a is a positive squeezing matrix and 
C = (fi,/7r) da / 2 Vdet A a is normalization constant. The 
purity can now be written as an integral 



I{t) « h~ M Jdjdj exp I lp(j) p (j), 
* = v(j A , J*b) - «(ja, Jb) + vQ A , Jb) - w(j A , Jb)- (8) 



St 



Next we expand the phase $ around the position j* = 
Ua.' Jb) °f ^ ne initial packet. The constant and the linear 
terms cancel exactly and the lowest order non- vanishing 
term is quadratic 



*«(jA-jA)-w AB (j*)(jB-j B ) 



(9) 



where v AB is a g?a x ^b matrix of mixed second derivatives 
of v evaluated at the position of the initial packet, 



( v AB)kl 



d 2 v 



9(jA)fc9(j B )i 



(10) 



Using this expansion in the integral for purity we see 
that the resulting 2d dimensional integral is Gaussian 
and can therefore be expressed in terms of a determinant 
of a 2d x 2d matrix. Using special properties of the re- 
sulting matrix the determinant can be reduced |Tlj to a 
determinant of a c?a x c?a matrix, with the final result 



I(t) 



v/det(l + (^) 2 u) 



« = a a 1s abA 



B V BA 



, (11) 



where u is a c?a x ^a matrix involving v AB and its trans- 
pose Vg A . Note that the matrix u is a classical quantity 
(independent of h) that depends only on observable v 
and on the position of the initial packet. This explicit 
formula for purity decay is the main result of the present 
paper jl5l |. 

Before discussing its consequences let us remind on its 
range of validity. The restrictions are rather weak: v 
must be nonvanishing (typical for regular systems) and 
smooth on a scale of the initial packet oc y%, time must 
be larger than the averaging time t > i avc , and the 
coupling must be small S < 1. In addition, the phase 
<!> should increment by a small amount for neighboring 
quantum numbers, which translates into the condition 

<^IKbII < 

The most prominent feature of the formula pi(l for 
the purity decay for initial product wave packets is 
its h independence. In the linear response calculation 
this ^-independence has already been theoretically pre- 
dicted 12| as well as numerically confirmed [13|. Here 
we have a full expression to all orders. We also see that 
the scaling of the decay time td on which I(t) decays is 
td ~ 1/8. This means that the purity will decay on a 
very long time scale and so the wave packets are univer- 
sal pointer states i.e. the most robust states. For 
small St we can expand the determinant and we get ini- 
tial quadratic decay I(t) = 1 — ^(St) 2 tr[u] + • • • . For 
large times we use the fact that det(l + zu) is a polyno- 
mial in z of order r = rank (u) , so we have asymptotic 
power law decay I(t) X const (8t)~ r . Note that the rank 
of u is bounded by the minimal of the subspace dimen- 
sions, i.e. 1 < r < min{dA,dBj-, since the definition © 
is symmetric with respect to interchanging the roles of 
the subspaces "A" and "B". Let us give two simple ex- 
amples: (i) For dA — 1 and for any g?b we will always 
have asymptotic power law decay with r = 1. If a single 
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DOF of the subsystem "A" is coupled with all DOF of 
the subsystem "B" , e.g. v = j A ® {j B \ + j'b2 + •■•)> trien 
\v"\ 2 oc d B and we have I(t) x l^St^/d^); (ii) Let us 
consider a multidimensional system where the matrix u 
is of rank one so it can be written as a direct product of 
two vectors, u = x (g) y. The determinant occurring in 



I(t) is then simply det (1 + (St) 2 u) = 1 + (5t) 2 x ■ y. Such 
is the case for instance if we have a coupling of the same 
strength between all pairs of DOF. The dot product is in 
this case x • y oc d A d B and we have I(t) X l/(Sty/d A d B ), 
i.e. the power of the algebraic decay is independent of 
both d\ and d B . 

In Q the author predicted a universal asymptotic t~ 2 
decay of purity independent of the dimensions or the cou- 
pling involved. He obtained this result for an average pu- 
rity, i.e. averaged over the position of the initial packet. 
Our result for I(t) clearly can not reproduce the re- 
sult of Ref. 0, even if we average over the position of 
the initial packet. Such an average decay will in general 
depend on the functional dependence of the matrix u on 
the position of the initial packet. 

We continue with a numerical demonstration of the 
theoretical prediction for purity decay ■ For the first 
example we take a 1 + 1 DOF system, d A = d B = 1, of two 
an-harmonic oscillators with the uncoupled Hamiltonian 



H a = lA {ha\a A - A) 2 + lB (ha+a B - A) 2 , (12) 

where a + ,a are standard boson raising/lowering opera- 
tors. For the coupling we take 

V = h 2 (a+ + a A ) 2 (a+ + a B ) 2 . (13) 

The corresponding classical Hamiltonian h reads 

h = 7a0a - A) 2 + 7 bOb - A) 2 + WSj A j B sin 2 6 A sin 2 6 B 

where 9 a are the canonical angles. The initial wave packet 
on both subsystems is a boson coherent state 



|Va(0)> = |Vb(0)) = \a) = e° 



*°|0), (15) 



where |0) is the ground state. The parameter a is chosen 
as a = \J j* /h with j* = 0.1. The squeezing parameter 
for the coherent states lfT5|l is Aa,b = l/(2j*). Other 
parameters of the Hamiltonian are 7a = 1, 7b = 0.6456. 
The offset A — 1.2 was chosen in order to have nonzero 
classical frequency dh/dj at the position of the initial 
packet. This is needed in order for v to be well defined. 
Time averaged coupling is calculated easily, v — 4jaJb- 
The matrix u is now just a number, u = (8j*) 2 . Theo- 
retical prediction for the purity decay is thus 



lit) 



1 



(16) 
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FIG. 1: Purity decay for 1+1 DOF system (11211^ for 6 = 0.04 
and different l/H = 10,25,50,100,500, from top to bottom. 
Dashed line with pluses is the theoretical formula 




The results of numerical simulation together with the 
theory are shown in Figs. ^ and [21 In Fig. ^ we see 



FIG. 2: Purity decay for 1 + 1 DOF system I|12ll3fl for h = 
1/100 and different 5 = 0.64, 0.32, 0.16, 0.08, 0.04, from left to 
right. Dashed lines give the theoretical prediction 116^ . 



that the decay is indeed ^.-independent, apart from a fi- 
nite size fluctuating plateau after long time. The size of 
this plateau is of the order I(t — > oo) ~ 1/N c s, where 
NeB ~ \/8j*/fi is an effective Hilbert space dimension, 
i.e. the number of action eigenstates overlaping with the 
initial coherent state l|15|) • Strong revivals for large h are 
a consequence of small number of available states N e g 
and low dimensionality. Revivals are expected to be less 
pronounced for larger dimensionalities d A , d B , similarly 
as for the fidelity For large times one can clearly 

observe asymptotic i _1 decay of the purity. In Fig. [3 we 
fix h and change the coupling strength 8 instead. Apart 
from oscillations we see a good agreement with the the- 
ory also for large 5. Oscillations for times t < 10 are a 
consequence of the fact that the time averaging of V © 
converges only after some averaging time i avc which is of 
order ~ 10 in our case. 

As for the second numerical example we take a 2 + 2 
DOF system (g?a.b = 2) which is the simplest case where 
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0.01 




500 



FIG. 3: Purity decay for a 2 + 2 DOF system and two 
different couplings showing different asymptotic power law 
decay. Full line is the numerics while the two dashed lines 
with crosses are theoretical predictions, the one with a smaller 
slope for (1191 and the other for 1181 . 



we can find different power of the asymptotic decay, de- 
pending on the topology of the coupling. The uncoupled 
Hamiltonian now reads 



H = -fi(hafax 



A)(fi-a^"a2 



A) 
A). 



(17) 



Subscripts "1" and "2" describe two DOF of the sub- 
system "A", while "3" and "4" compose the subsys- 
tem "B". The parameters are 71 = 1, 72 = 0.64 and 
A = 1.2. The initial state is a product state of four 
boson coherent states, |^i, 2,3,4(0)) = \ct), all with the 
same a = \J j* /fi. For the coupling we consider two 
cases which will give different power of the asymptotic 
decay. Case I.: V = V13 + V24, where two coupling terms 
are of the same form as for 1 + 1 DOF system l|13f) and 
the indices denote between which two degrees of freedom 



the coupling acts. The matrix u as well as the relevant 
determinant is easily calculated resulting in a simple ex- 
pression for the purity l)ll|l 



I(t) = T 



1 



(8j*sty 



(18) 



We see that we have a quadratic asymptotic decay, I{t) >c 
l/(St) 2 . Case II.: All to all coupling, V = V13 + Vu + 
V23 + V24, results in a rank-one (r = 1) matrix u giving 
the purity decay lfTT|) 



/(*) = 



1 



y/l + (16j*(Ji) 2 



(19) 



Results of the numerical simulation for both cases are 
shown in Fig. [3] The coupling strength and the location 
of the initial packets are 8 = 0.04, j* — 0.1 for the case 
I., and S = 0.02, j* = 0.2 for the case II. From Fig. 
we see that one indeed has asymptotic t^ 1 or i -2 decay, 
depending on the topology of the coupling. 

In conclusion, we have derived purity decay for ini- 
tial localized wave-packets in bipartite systems with a 
non-vanishing (non-trivial) time-averaged coupling oper- 
ator. Such situation naturally occurs in systems where 
an uncoupled part of the Hamiltonian represents regular 
dynamics. Purity decays in time inversely proportional 
to the coupling strength and is independent of Planck's 
constant. The decay is algebraic with the asymptotic 
power-law exponent ranging between 1 and the minimal 
dimension of the subsystems depending on the topology 
of the coupling. 
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